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基于 人 神经 网 络 的 半 接 双星 测 光 解 轨 方法 
REF, aa, AR, FAH, aka? 
(1. 中 国 科学 院 国 家 天 文 台 北京 100101; 2. 中 国 科学 院 大 学 北京 100049) 


摘 要 : 半 接 双星 对 于 研究 相互 作用 双星 的 形成 和 演化 有 着 重要 意义 。 随 着 大 规模 时 域 巡 天 
时 代 的 到 来 ， 预 计 将 发 现 大 量 此 类 天 体 。 针 对 海量 的 时 域 观测 数据 ， 需 要 一 种 快捷 的 建 模 工具 
对 半 接 双星 开展 自动 化 光 变 曲线 分 析 。 基 于 神经 网 络 构建 了 半 接 双星 快速 测 光 解 轨 模 型 。 该 模 
型 根据 光 变 曲线 和 已 知 的 主星 温度 对 半 接 双星 的 轨道 进行 建 模 ， 获 得 其 轨道 倾角 、 相 对 半径 、 
质量 比 和 温度 比 等 4 个 基本 参数 。 结 果 表 明 ， 半 接 双 星 神经 网 络 解 轨 模 型 能 够 快速 地 对 一 条 光 变 
曲线 建 模 。 在 测 光 误差 小 于 光 变 曲线 振幅 1% 时 ， 模 型 对 于 轨道 倾角 接近 90" 、 温 度 比 约 为 0.6、 


光 变 振幅 为 1.84mag 的 半 接 双星 的 轨道 倾角 、 相 对 半径 、 质量 比 和 温度 比 的 测量 误差 分 别 为 1.251、 


0.004、0.008 和 0.003。 另 外 ， 模 型 应 用 在 Kepler 实 测 光 变 曲 线 的 结果 表明 ， 模 型 能 够 较为 准确 地 
对 脉动 食 双星 的 光 变 曲线 进行 建 模 ( 拟 合 度 可 达 0.9 以 上 ) 。 此 外 ， 该 模型 作为 一 个 通用 工具 可 
以 迁移 到 不 同 测 光 巡天 项 目 上 。 
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双星 系统 为 研究 恒星 结构 与 演化 提供 了 一 个 重要 信息 来 源 ， 双 星 中 的 半 接 双星 (semi- 
detached binary)， 是 指 双星 中 的 一 颗 子 星 充 满 洛 希 办 ， 另 外 的 一 颗 子 星 未 充满 洛 希 办 ， 并 且 充 
满 洛 希 办 的 子 星 通过 内 拉 格 朗 日 点 向 另外 一 颗 子 星 传递 物质 的 双星 系统 。 在 双星 系统 中 ， 子 星 
充满 洛 希 办 发 生物 质 转移 ， 会 引起 双星 系统 的 动能 、 势 能 、 热 能 的 变化 ， 这 将 直接 影响 双星 的 
结构 和 演化 ， 同 时 这 会 使 得 恒星 演化 进程 和 结局 发 生 戏 剧 性 的 改变 趾 。 有 关 相 互 作用 双星 的 演 
化 仍然 是 目前 双星 演化 理论 中 的 难点 外 。 半 接 双 星子 星 间 的 相互 作用 和 两 子 星 间 的 物质 交流 和 
能 量 转移 等 过 程 ， 使 得 半 接 双星 的 演化 进程 和 结局 变 得 难以 确定 ， 因 此 ， 相 互 作 用 的 双星 的 形 
成 、 演 化 和 结局 仍 是 双星 演化 的 未 解 问 题 丫 。 半 接 双 星 的 参数 (如 质量 比 、 轨 道 倾角 和 温度 比 
等 ) 为 研究 密 近 双星 系统 的 质量 和 角 动 量 转移 、 损 失 和 重新 分 布 等 提供 了 有 利 的 数据 基础 。 除 
此 之 外 ， 光 变 昌 线 的 建 模 也 是 研究 食 双星 中 恒星 脉动 机 制 的 基础 ， 通 过 光 变 曲线 建 模 消 除 双星 
掩 食 造 成 的 光度 变化 后 ， 我 们 能 够 发 现 双 星 中 的 非 径 向 振荡 ， 进 而 通过 星 震 分 析 来 约束 超 射 参 
数 ， 为 了 解 恒星 内 部 结构 提供 重要 依据 "1。 

近年 来 ， 随 着 NASA/Kepler 和 TESSW 空 间 观测 站 的 发 射 ， 以 及 ZTFÜ'RI ASAS-SNDI 等 望 
远 镜 的 成 功 运 行 ， 我们 进入 了 一 个 精确 空间 大 规模 时 域 测 光 的 新 时 代 ， 越 来 越 多 的 双星 也 随 之 
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被 证 认 ，Kepler/K2 测 光 巡天 证 认 了 约 2900 颗 食 双星 0 中, 文 [11] 从 ZTF DR2 数据 集中 证 认 并 提 
供 了 约 350,000 颗 食 双 星 的 光 变 曲线 ,ASAS-SN 测 光 巡 天 发 布 了 约 426,000 颗 变星 的 光 变 曲线 
Catalina 项 目 提 供 了 约 47,000 条 周期 变星 的 光 变 曲线 小。 这 为 双星 的 发 现 和 研究 提供 了 巨大 的 
机 遇 。 目 前 ， 对 于 双星 的 光 变 曲线 建 模 大 多 基于 PHOEBE 模型 或 是 Wilson-Devinney(WD) 
codel", 采用 WD 或 者 PHOEBE 程序 对 一 个 目标 的 光 变 曲线 建 模 需 要 的 计算 时 间 长 达 数 小 时 其 
至 数 天 。 海 量 的 光 变 曲线 需要 一 种 更 加 快速 的 光 变 曲线 建 模 方法 ， 这 对 开展 双星 的 性 质 统计 分 
析 以 及 快速 搜寻 双星 中 具有 特殊 机 制 的 样本 有 重要 意义 。 

本 文 提 出 采用 神经 网 络 对 半 接 双星 的 光 变 曲线 进行 轨道 参数 求解 ， 达 到 光 变 曲线 的 快速 建 
模 的 目的 。 


1 测 光 解 轨 的 神经 网 络 建 模 


一 般 来 说 ， 双 星 的 光 变 曲线 中 食 的 形状 可 以 反映 轨道 倾角 ， 掩 食 的 持续 时 间 和 它们 的 相对 
S 位 置 反映 了 轨道 的 偏心 率 、 近 星 点 经 度 和 两 个 子 星 的 半径 之 和 ， 除 此 之 外 ， 两 次 掩 食 的 深度 比 
LO 表明 了 两 颖 子 星 表面 温度 和 遮挡 面积 之 间 的 关系 ， 另 外 ， 重 力 昏 暗 、 临 边 变 暗 、 反 射 系数 、 黑 
li 子 以 及 恒星 周围 的 物质 也 会 在 光 变 曲线 上 有 一 定 的 显示 。 由 此 ， 通 过 解 轨 ， 我 们 可 以 从 光 变 曲 
= 线 中 得 到 双星 系统 的 周期 、 轨 道 倾角 、 光 度 比 、 相 对 半径 Ra, 其 中 R 为 子 星 的 等 效 半径 ， 
一 a 为 双星 轨道 半 长 轴 ) 等 特征 。 近 相 接 双星 光 变 曲线 一 些 特征 的 统计 显示 ， 轨 道 倾 角 大 于 70° 的 
近 相 接 双星 才 容 易 观测 ， 并 且 大 部 分 近 相 接 双星 是 由 光谱 型 为 A，F 型 的 主星 和 晚 一 到 两 个 光 
谱 型 的 次 星 组 成 [9 。 
c 1. 1 模型 样本 数据 
N 半 接 双星 根据 子 星 是 否 充满 洛 希 准 ， 半 接 双 星 充满 络 希 准 的 类 型 可 以 分 为 两 种 09，(1) E 
2 星 (质量 大 的 子 星 ) FSA, ULATION RIES AE, DKE EMITE) 充 
X HARI TELA TORII IRE. 
c 由 于 我 们 无 法 直接 从 光 变 曲线 判定 半 接 双星 是 主星 充满 还 是 次 星 充满 的 状态 ， 并 且 两 种 状 
rz 态 下 的 光 变 曲线 在 某 些 参数 空间 存在 简 并 。 如 果 将 主星 充满 和 次 星 充满 的 光 变 曲线 进行 统一 训 
O 练 ， 由 于 参数 简 并 导致 模型 精度 降低 ， 造 成 模型 预测 的 参数 不 准确 ， 且 无 法 判定 半 接 双星 的 状 
态 。 因 此 ， 本 文 分 别 构建 主星 充满 和 次 星 充满 两 个 状态 的 光 变 曲线 测 光 解 轨 模 型 。 针 对 一 条 实 


测 光 变 曲线 ， 模 型 可 以 给 出 在 主星 充满 和 次 星 充满 两 个 状态 下 的 最 佳 参数 解 ， 并 由 PHOEBE 根 
据 两 组 参数 重 现 的 光 变 曲线 联合 实测 光 变 曲 线 计算 拟 合 度 ， 最 后 根据 拟 合 度 判 断 半 接 双星 所 处 
的 状态 以 及 参数 。 

我 们 将 主星 温度 T 的 范围 设 定 为 5000 - 10000K; 并 将 “主星 ”定义 为 双星 中 质量 较 大 的 子 
星 ， 并 且 “ 主 星 ” 也 定义 为 “ 星 1”。 轨 道 倾角 (incD) 的 范围 为 60 - 90"， 质 量 比 (q=) 的 范围 为 
0.1 - 1， 次 星 温度 (7?) 的 范围 假设 为 3800 - 10000K， 子 星 洛 希 办 的 临界 半径 R11 和 R11 "为 
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并 且 在 该 范围 内 随机 生成 相对 半径 (R/a)。 
基于 上 述 的 参数 范围 ， 和 针对 主星 充满 测 光 解 轨 模型 ， 我 们 将 轨道 倾角 、 质 量 比 、 主 次 星 温 
度 、 次 星相 对 半径 作为 自由 参数 ， 在 规定 的 各 个 参数 范围 内 每 次 均匀 、 随 机 地 抽取 一 组 参数 ， 
利用 PHOEBE 在 0-1 的 相位 内 ， 生 成 100 个 采样 点 的 光 变 曲线 ， 并 将 光 变 曲线 的 流量 值 转换 成 
Kepler 星 等 ， 从 而 获得 构建 主星 充满 测 光 解 轨 模型 的 数据 集 。PHOEBE 用 “流量 (flux)” 表 示 双 
星 的 光 变 曲线 的 振幅 ， 因 此 需要 根据 
mag; = —2.5log10(flux;) ， (3) 
mag, = magi 一 Zak (4) 
对 PHOEBE ERICA HRTEM. Hon. flux HAHH PHOEBE 生成 100 个 采样 点 的 光 
变 曲 线 的 流量 值 ; mag, 为 标准 化 后 的 光 变 曲线 对 应 的 星 等 。 同 样 ， 针 对 次 星 充 满 模型 采取 同 
样 的 方法 构建 模型 的 数据 集 。 
针对 主星 充满 测 光 解 轨 模型 和 次 星 充满 测 光 解 轨 模 型 , 本 文 分 别 产生 了 约 30 万 条 光 变 曲线 
用 于 构建 模型 。 图 1(a) 和 (b) 分 别 为 数据 集 在 各 个 参数 的 频率 统计 直方 图 ， 从 图 1 可 以 看 出 ， 数 
据 集 在 相对 半径 、 温 度 比 上 的 分 布 较 不 均匀 ， 这 是 由 于 参数 之 间 的 相互 约束 导致 的 。 因 此 ,为 
了 解决 样本 在 各 个 参数 空间 分 布 不 均匀 的 问题 ， 对 于 模型 训练 中 每 次 迭代 过 程 随机 抽取 训练 集 
中 的 5 万 条 数据 进行 训练 ， 直 至 模型 收 伍 。 同 时 ， 数 据 集 中 的 主星 温度 、 轨 道 倾角 、 质 量 比 、 
主 次 星 温度 比 、 次 星相 对 半径 的 尺度 也 通过 skleam 提供 的 StandardScaler 标准 化 方法 进行 了 标 
准 化 。 最 终 , 将 数据 集 按照 训练 集 : 验证 集 : 测试 集 =7: 2: 1 进行 划分 用 于 模型 的 训练 和 检验 。 
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图 1 半 接 双星 测 光 解 轨 模 型 的 数据 集 在 各 个 参数 上 的 频率 统计 直方 图 。(a) 主 星 充满 b) 次 星 充满 


Fig. 1 The frequency distribution of the five parameters (T1, incl, q, R/a, T2/T1) for light curve modeling of semi-detached binaries. (a) 


The primary filling model; (b) the secondary fill model. 


1.2 测 光 解 轨 的 神经 网 络 建 模 
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图 2 显示 了 本 文采 用 的 网 络 结构 。 本 文 基于 多 层 感 知 器 (Multi-Layer Perception, MLP) 结 
构 设 计 模 型 ， 多 层 感知 器 包括 输入 层 、 隐 藏 层 和 输出 层 。 本 文 的 网 络 结构 包含 1 个 输入 层 ，8 个 
隐藏 层 和 1 个 输出 层 。 输 入 层 包 含 101 个 神经 元 , 即 100 个 数据 点 构成 的 光 变 曲线 和 主星 温度 (7T1)。 
每 个 隐藏 层 是 一 个 全 联接 层 (FC，full connection)， 每 层 输入 和 输出 的 神经 元 个 数 如 图 2 中 的 数字 
所 示 。 为 了 减 小 传播 过 程 中 的 特征 损失 ， 在 隐藏 层 ， 本 文通 过 引入 残 差 块 增强 特征 。 残 差 块 的 
隐藏 层 之 间 采 用 tanp 激 活 函数 进行 连接 ， 有 效 增强 了 网 络 非 线性 表达 能 力 。 本 文采 用 反 向 传播 
算法 对 多 层 感 知 模型 进行 训练 ， 采 用 Adam 优 化 器 和 损失 函数 L2， 即 

Dia = FL: — f)? (5) 
其 中 ，y; 为 真实 值 ，f (xi) 为 预测 值 。 最 后 ， 模 型 输出 轨道 倾角 (incl)、 相 对 半径 (R/a)、 质 量 比 (g) 
和 温度 比 (TT) 4 个 参数 。 
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图 2 网 络 结构 图 


Fig.2 Construction of the neural network for light curve modeling 


1.2.1 构建 主星 充满 测 光 解 轨 模 型 


在 构建 主星 充满 测 光 解 轨 模 型 时 , 我 们 将 训练 集中 100 个 采样 点 构成 的 光 变 曲线 和 主星 温度 
(7T1) 作 为 输入 数据 , 将 轨道 倾角 (incl)、 次 星相 对 半径 (Ry/a)、 质量 比 (q)、 温度 比 (72/7T1) 作 为 输出 ， 
对 光 变 和 参数 之 间 的 映射 关系 进行 了 训练 ,结果 如 图 3, 图 3(a)-(d) 依 次 显示 了 测试 集 的 轨道 倾角 、 
次 星相 对 半径 、 质 量 比 、 温度 比 和 模型 预测 得 到 的 结果 对 比 。 从 图 3 可 以 看 出 模型 预测 的 轨道 倾 
角 、 次 星相 对 半径 、 质 量 比 温度 比 和 测试 集 的 真实 值 的 偏差 分 别 为 0.039, 0.000, -0.003 和 -0.002， 
并 且 轨 道 倾角 、 相 对 半径 、 质 量 比 和 温度 比 与 测试 集 的 真实 值 的 残 差 的 标准 差分 别 为 0.238, 0.003, 
0.010 和 0.024。 
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图 3 主星 充满 的 测 光 解 轨 模 型 进行 检验 的 结果 。(a): 测 试 集 的 轨道 倾角 (ncl) 和 模型 预测 得 到 轨道 倾角 的 结果 对 比 ;(b): 测 试 集 的 次 星相 
对 半径 (Rya) 和 模型 预测 得 到 相对 半径 的 结果 对 比 ; (ec) 测试 集 的 质量 比 (q) 和 模型 预测 得 到 质量 比 的 结果 对 比 ; (qd) 测 试 集 的 温度 比 (Ty/T) 
和 模型 预测 得 到 的 温度 比 结果 对 比 。 
Fig. 3 The result of the primary filling model. (a)-(d) show the residual between model prediction and the true value of the test set in orbital 
inclination(incl), relative radius(R/a), mass ratio(q), and temperature ratio(T2/T1). 


图 4 显示 了 假设 主星 充满 并 根据 主星 温度 TI=6000K， 轨 道 倾角 incl=85°*， 次 星相 对 半径 
Rs/a=0.3， 质 量 比 g=0.5， 温 度 比 TyT1=0.866 为 参数 模拟 的 一 条 半 接 双星 光 变 曲线 利用 主星 充满 
测 光 解 轨 模型 进行 光 变 曲线 建 模 的 结果 。 模 型 根据 输入 的 光 变 曲线 和 主星 温度 ， 预 测 得 到 的 轨 
道 倾 角 为 84.9$"， 相 对 半径 为 0.299， 质 量 比 为 0.499， 温 度 比 为 0.867。 接 着 ， 我 们 根据 模型 预测 
得 到 的 参数 利用 PHOEBE 对 光 变 曲线 进行 了 重 现 ， 结 果 如 图 4 上 方 图 像 中 的 红色 实 线 ， 图 4 的 下 
方 图 像 显示 了 原始 光 变 曲线 《〈 蓝 色 数 据点 ) 和 采用 模型 预测 得 到 的 参数 重 现 的 光 变 曲线 之 间 的 
残 差 。 从 图 中 可 以 看 出 ， 两 条 光 变 曲线 能 够 较 好 地 匹配 ， 而 在 掩 食 对 应 位 置 光 变 曲线 的 残 差 大 
于 其 他 位 置 的 残 差 , 主要 是 模型 在 预测 时 输入 的 100 个 数据 点 的 光 变 曲线 是 均匀 采样 ,使 得 掩 食 
部 分 的 数据 点 较 少 造成 的 。 此 外 ， 模 型 预测 得 到 的 参数 重 现 的 光 变 曲 线 和 原始 观测 数据 之 间 的 
残 差 小 于 0.005mag， 并 且 拟 合 度 (R?) 为 0.994。 其 中 拟 合 度 定义 为 


n S-yy? V 3 mx. 
R?=1- Ae (f(xi) 为 预测 值 ，y 为 真实 值 ) . (6) 
i=1W1 


* tnetc —— DL pre 
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图 4 采用 主星 充满 测 光 解 轨 模型 进行 光 变 曲线 建 模 的 结果 ,其 主要 参数 为 主星 温度 TI=6000K， 轨 道 倾角 incl=85°*， 相 对 半径 Rya=0.3， 
质量 比 g=0.5， 温 度 比 7y7T1=0.866 
Fig.4 The validation of the light curve for the primary filling model, the blue dots are the true light curve, and the red line is the model 
prediction result. The parameters are TI=6000K, incl=85°, Rz/a=0.3, q=0.5, and T/T1=0.866. 
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1. 2.2 构建 次 星 充满 测 光 解 轨 模型 

同样 的 ， 我 们 将 训练 集中 100 个 采样 点 构成 的 光 变 曲线 和 主星 温度 (7 作为 输入 数据 ， 将 轨 
道 倾角 、 主 星相 对 半径 、 质 量 比 和 温度 比 作为 输出 ， 对 次 星 充 满 测 光 解 轨 模型 中 的 光 变 曲线 和 
参数 之 间 的 映射 关系 进行 了 训练 。 结 果 如 图 5， 图 5(a)-(qd) 依 次 显示 了 测试 集中 轨道 倾角 、 相 对 半 
径 、 质 量 比 、 温 度 比 和 模型 预测 得 到 的 结果 对 比 ， 从 图 中 可 以 看 出 ， 模 型 预测 的 轨道 倾角 、 相 
对 半径 、 质 量 比 、 温 度 比 和 测试 集 的 真实 值 的 偏差 分 别 为 0.013，0.000，0.002 和 -0.001， 并 且 轨 
道 倾角 、 相 对 半径 、 质 量 比 和 温度 比 与 测试 集 真实 值 残 差 的 标准 差分 别 为 0.291，0.004，0.013 
和 0.016。 


man= 20 wtd-0 004 


图 5 次 星 充 满 的 测 光 解 轨 模 型 进行 检验 的 结果 。(a): 测 试 集 的 轨道 倾角 (ncD 和 模型 预测 得 到 轨道 倾角 的 结果 对 比 ; (b): 测 试 集 的 主星 相 
对 半径 (Rya) 和 模型 预测 得 到 相对 半径 的 结果 对 比 ; (c) 测 试 集 的 质量 比 (q) 和 模型 预测 得 到 质量 比 的 结果 对 比 ; (d) 测 试 集 的 温度 比 (Tw/T1) 
和 模型 预测 得 到 的 温度 比 结果 对 比 。 
Fig.5 The result of the secondary filling model. (a)-(d) show the residual between model prediction and the true value of the test set in 
orbital inclination(incl), relative radius(R/a), mass ratio(q), and temperature ratio(T2/T1) 


同样 ， 假 设 次 星 充满 ， 图 6 显示 了 主星 温度 娓 =7600K， 轨 道 倾 角 imci=85"， 主 星相 对 半径 
Ri/a-0.4, 质量 比 q=0.2, 温度 比 Tz/T=0.684 的 一 条 光 变 曲线 采用 次 星 充满 模型 进行 建 模 的 结果 。 
模型 根据 输入 的 光 变 曲线 和 主星 温度 ， 获 得 的 参数 为 轨道 倾角 为 84.88"， 相 对 半径 为 0.399， 质 
量 比 为 0.200， 温 度 比 为 0.681。 图 中 蓝 色 点 显示 了 原始 模拟 的 光 变 曲线 ， 红 色 实 线 是 PHOEBE 根 
据 模 型 预测 的 参数 重 现 的 光 变 曲线 ， 从 图 6 可 以 看 出 ， 两 条 光 变 曲线 能 够 较 好 地 匹配 ， 与 主星 充 
满 测 光 解 轨 模 型 一 样 ， 图 6 中 的 光 变 曲 线 在 掩 食 对 应 位 置 的 残 差 也 大 于 其 他 位 置 的 残 差 , 这 主要 
也 是 因为 掩 食 部 分 的 数据 点 较 少 ,另外 ,从 图 中 可 以 看 出 两 条 光 变 曲线 之 间 的 残 差 小 于 0.005mag， 
并 且 拟 合 度 (R?) 为 0.9957。 


» 
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图 6 采用 次 星 充 满 测 光 解 轨 模 型 进行 光 变 曲线 建 模 的 结果 ， 其 主要 参数 为 TI=7600K，incl=85°*，Ri/a=0.4， 质 量 比 g=0.2， 温 度 比 
T/T1=0.684 


Fig.6 The validation of the light curve for the secondary filling model, the blue dots are the true light curve, and the red line is the model 
prediction result. The parameters are TI=7600K, incl=85°, Ri/a=0.4, q=0.2, and T2/T)=0.684. 


3 实测 数据 的 测 光 解 轨 分 析 


3.1 Kepler 天 区 半 接 双星 测 光 解 轨 结 果 

NASA 的 Kepler 卫 星 于 2009 年 3 月 发 射 ， 其 主 镜 口径 为 1.4m( 有 效 通 光 口径 为 0.9$m)， 视 场 为 
105 平 方 度 。Kepler 卫 星 提供 了 长 上 曝光 Clongeadence, LC; 曝光 时 间 29.4min) 的 观测 数据 和 短 
曝光 (shortcadence，SC; 曝光 时 间 58.8s) 的 观测 数据 。 文 [18] 对 Kepler 卫 星 观 测 得 到 的 2165 颗 食 
双星 进行 分 类 ， 共 有 1261 个 分 离 双星 (detached binary)，152 个 半 接 双星 (semi-detached binary)， 
469 个 相 接 双 星 (overcontact binary)。 我 们 从 其 中 挑选 已 经 获得 轨道 参数 解 的 半 接 双星 对 模型 应 用 
在 实际 观测 的 数据 上 的 效果 进行 了 分 析 。 

针对 实际 的 观测 数据 ， 解 轨 步 骤 如 下 : 

C1) 光 变 曲线 预 处 理 ， 首 先 ， 将 光 变 曲线 根据 Kepler 卫 星 提供 的 周期 折 县 到 相位 上 ， 然 后 
根据 (2) 式 和 (3) 式 对 观测 数据 进行 标准 化 并 剔除 观测 误差 大 于 3omag 的 数据 点 ; 接着 ， 我 们 将 原 
始 数 据 在 相位 区 间 (0,1) 上 分 为 100 个 bin， 并 将 每 个 bin 中 数据 点 的 均值 作为 输入 的 光 变 曲线 。 图 
7(a)-(e) 显 示 了 Kepler 卫星 观测 获得 的 5 颗 半 接 双 星 CKICIOS81918"!, KIC10619109"! , 
KIC06669809"!, KIC06206751!'FIKIC060481067") 的 光 变 曲线 经 过 预 处 理 后 的 结果 ， 其 中 灰 
色 点 表示 原始 的 观测 数据 ， 蓝 色 数 据点 是 经 预 处 理 后 输入 模型 中 进行 参数 预测 的 光 变 曲线 。 
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图 7 KIC10581918, KIC10619109, Clee KIC06206751 和 KIC06048106 经 过 预 处 理 步 又 的 结果 。 灰 色 数 据点 是 观测 获 
得 的 光 变 曲线 ， 蓝 色 数 据点 是 经 过 预 处 理 后 输入 到 横 型 中 的 光 变 曲线 

Fig.7 The results of light curve pre-processing for KIC10581918, KIC10619109, KIC06669809, KIC06206751, and KIC06048106. 

The Grey points are the light curves from observation, and blue dots are the result of pre-processing. 


(2) 将 第 1 步 获 得 的 光 变 曲线 以 及 由 Kepler 卫 星 提供 的 (或 从 光谱 测量 得 到 的 ) 主星 温度 娓 分 别 
输入 主星 充满 和 次 星 充满 测 光 解 轨 模 型 中 对 参数 预测 ， 并 将 预测 的 参数 输入 PHOEBE 中 进行 光 
变 曲 线 重 现 ， 计 算 拟 合 的 光 变 曲线 与 观测 数据 的 拟 合 度 RR， 最 后 将 RR 大 的 光 变 曲线 对 应 的 参数 
作为 测 光 解 轨 的 参数 。 

根据 上 述 步 又 ， 我 们 对 5 颗 半 接 双 星 进行 参数 预测 。 图 8(a)-(e) 分 别 为 采用 本 文 的 模型 对 
KIC10581918, KIC10619109, KIC06669809, n e. 

WAR, Bsp tU ex JUGAR ZR, RECS ER AJ KH] FE A Te DICERES BS s HR 
色 线 为 采用 次 星 测 光 解 轨 模 型 得 到 的 结果 。 从 图 中 可 以 看 出 ,5 颗 半 接 双星 都 是 采用 次 星 充满 模 
型 预测 得 到 的 参数 拟 合 的 光 变 曲线 与 实际 观测 的 光 变 曲线 拟 合 度 高 于 主星 充满 的 模型 预测 得 到 
的 结果 ,因此 ,我 们 认为 KIC10581918, KIC10619109, KIC06669809, KIC06206751 和 KIC06048106 
是 次 星 充满 的 半 接 双星 ， 这 与 已 有 文献 得 到 的 结果 一 致 。 
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图 8 (a)-(e) 依 次 显示 模型 拟 合 KIC10581918，KIC10619109，KIC06669809，KIC06206751 和 KIC06048106 的 结果 ， 灰 色 数 据点 表 
示 Kepler 观 测 获 得 的 光 变 曲线 ， 橙 色 虚 线 显示 的 是 采用 主星 充满 测 光 模型 求解 的 参数 拟 合 得 到 的 光 变 曲线 ， 蓝 色 线 显示 的 是 采用 次 星 
测 光 解 轨 模 型 求解 的 参数 拟 合 得 到 的 光 变 曲线 
Fig.8 The results of model prediction for KIC10581918, KIC10619109, KIC06669809, KIC06206751, and KIC06048106. The Grey 
points are the light curves from observation, the blue lines are the result of the secondary filling model, and the orange dashed dot lines result 
from the primary filling model. 


表 1 显 示 了 采用 基于 神经 网 络 的 半 接 双星 测 光 解 轨 模 型 获得 的 KIC10581918, KIC10619109, 
KIC06669809，KIC06206751 和 KIC06048106 的 参数 与 己 有 文献 中 采用 PHOEBE 或 WD 获 得 的 轨 
道 参数 的 对 比 。 从 表 1 可 以 看 出 , 采用 本 文 提 出 的 基于 神经 网 络 的 测 光 解 轨 模 型 对 半 接 双星 进行 
解 轨 得 到 的 参数 大 多 在 文献 中 获得 的 轨道 参数 的 误差 范围 内 , 对 于 KIC06669809 和 KIC06206751 


其 拟 合 度 低 于 0.9， 可 能 的 原 


在 0.25 和 0.75 附 近 ) 展示 不 对 称 的 
PK; (2) 对 于 KIC06669809， 除 了 交 


4 


变 


因 是 (1)KIC06669809 和 KIC06206751 的 光 变 


表 况 ， 这 可 能 是 由 于 
| 线 不 对 称 ， 


子 星 上 的 
其 轨道 倾角 也 低 于 7 


| 线 在 非 掩 时 刻 ( 相 位 


子 或 热 斑 造成 的 光 变 曲 线 不 对 


0”， 因 此 在 一 


定 程度 上 


影响 了 参数 的 估计 准确 度 ; (3) 对 于 KIC06206751， 除 了 在 相位 =0.25 和 0.75 附 近 光 变 曲 线 不 对 称 
的 原因 ， 还 有 可 能 是 由 双星 轨道 存在 一 定 的 偏心 率 导 致 的 (separation between primary and 


secondary eclipse (Seq)=0.5008)。 对 村 
于 没有 加 入 此 类 光 变 曲线 ， 因 此 模型 


于 没有 因子 星 上 的 黑子 、 热 斑 造 成 光 变 曲线 的 不 对 称 的 
这 类 食 双 星 的 光 变 曲线 进行 建 模 ， 这 对 后 续 分 析 恒 星 脉动 机 制 有 习 


该 模型 能 


子 星 上 有 黑子 或 热 斑 的 半 接 双星 ， 在 模型 构建 过 程 中 ， 由 
1 对 于 这 类 双星 的 光 变 曲线 建 模 ， 拟 合 度 可 能 低 于 0.9。 
脉动 食 双星 ， 该 模型 能 够 较为 准确 


EE 要 的 意义 。 男 外 ， 


而 对 
有 地 对 


够 在 对 大 量 光 变 曲线 进行 快速 建 模 后 ， 也 为 我 们 发 现 一 些 具有 特殊 机 制 的 食 双 星 提 供 了 一 个 有 
效 的 工具 。 
表 1 基于 深度 学 习 网 络 获得 的 参数 与 已 有 文献 结果 的 对 比 
Table 1 Comparison of the parameters between the deep learning results with the literature. 
Star Ti incl R/a q T/T R? 
Ref DL Ref DL Ref DL Ref DL 
KIC10581918" 7900 89.0+0.2 87.79 0.213+0.031 0.202 0.12+0.01 0.122 0.627+0.016 0.636 | 0.9701 
KIC10619109*! 7070 75.0+0.01 — 71.94 0.240+0.030 — 0.213 0.14+0.01 0.275 0.558+0.013 0.580 | 0.9348 
KIC06669809! 7400 67.5+0.1 71.70 0.399+0.046 — 0.417 0.310+0.010 0.164 0.661+0.018 0.696 | 0.8421 
K1C06206751?! 6965 75.000.200 | 7443 0.235+0.034 — 0.262 0.1295+0.0031 — 0.115 0.679+0.035 0.727 | 0.8979 
KIC06048106""! 6999 73.900.080 — 7228 0.218+0.068 — 0.199 0.2097+0.0033 — 0304 0.639+0.043 0.650 | 0.9238 
3. 2 模型 解 轨 结 果 分 析 
对 于 Kepler 巡 天 ， 针 对 Kpmag<l3mag 的 恒星 曝光 30min ， 其 测 光 精度 约 为 100ppm 
(~0.0001mag)” 1。 但 是 对 于 双星 ， 观 测 误差 (0 ，。 a) 对 食 双 星 的 参数 测量 的 影响 与 光 变 曲线 的 振 
幅 (4mag) 有 关 。 因 此 , 我 们 定义 相对 测量 误差 = o /Awag 来 表征 观测 误差 对 双星 的 光 变 曲 线 的 
影响 。 本 文 以 KIC10581918 为 样本 (振幅 约 为 1.84 mag) ， 利 用 蒙特 卡 洛 模拟 对 该 模型 在 不 同 观 
测 误差 下 参数 测量 的 误差 进行 估计 。 首 先 ， 我 们 对 KIC10581918 的 原始 观测 数据 加 上 不 同 相 对 


误差 对 应 的 噪声 ， 再 将 原始 数据 在 0-1 的 相 
为 输入 的 光 变 曲线 ， 利 用 次 星 充满 模型 估计 不 同 相 对 i 


对 半径 、 质 量 比 和 温度 比 以 及 它们 的 误差 。 图 
参数 与 已 有 文献 的 参数 的 差 值 ， 灰 色 阴 影 部 


时 , 模型 对 KIC10581918 的 轨道 倾角 、 相 对 半径 、 质 量 比 和 温 
0.008 和 0.003。 与 表 1 的 结果 相 比 ， 在 相对 误差 为 0.01 时 
度 比 的 测量 与 已 有 文献 的 结果 相 比 ， 
角 的 测量 ， 模 型 获得 的 轨道 倾角 
增 大 ， 模 型 预测 得 到 的 轨道 倾 


质量 比 和 温 


分 表示 模型 


浅 的 次 级 小 时 刻 的 影响 大 于 主 极 小 ， 


吴 差 下 KIC10581918 的 轨道 倾角 、 
9 中 横 坐 标 为 不 同 的 相对 误差 ,红色 点 为 模型 
J 预测 的 参数 误差 。 

从 图 9 可 以 看 出 ， 模 型 预测 参数 的 误差 随 着 相对 误差 的 增 大 而 增 大 ， 并 


Hii 


, 


度 比 的 测量 
针对 KIC10581918， 模 型 对 相对 半径 、 


与 已 有 文献 的 轨道 倾 
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图 9KIC10581918 在 不 同 观测 误差 下 ， 模 型 的 参数 测量 误差 的 统计 图 
Fig.9 The measurement error of orbital inclination(incl), relative radius(R/a), mass ratio(q), and temperature ratio(T2/T1) with 
different relative errors in observation for KIC10581918. 


综 上 可 得 ， 模 型 用 于 求解 Kepler 巡 天 的 半 接 双星 结果 表明 ， 该 模型 能 够 较为 准确 地 获得 双 
星 的 参数 , 并 且 根 据 模型 获得 的 参数 重 现 的 光 变 曲线 与 实测 的 光 变 曲线 的 拟 合 度 达到 0.9。 另外， 
在 测 光 误 差 相 对 于 光 变 曲线 的 振幅 为 0.01 时 ， 对 于 轨道 倾角 接近 90”、 温 度 比 约 为 0.6、 光 变 振 
幅 为 1.84mag 的 半 接 双星 ,根据 次 星 充满 测 光 解 轨 模型 获得 的 轨道 倾角 和 、 相 对 半径 、 质 量 比 和 温 
度 比 的 测量 误差 分 别 为 1.251，0.004，0.008 和 0.003， 这 为 批量 获得 半 接 双星 参数 研究 密 近 双星 
系统 的 质量 和 和 角 动量 转移 、 损 失 和 重新 分 布 等 提供 了 有 利 的 数据 基础 。 

该 神经 网 络 构建 的 半 接 双星 测 光 解 轨 模型 能 够 快速 地 对 一 条 光 变 曲线 建 模 , 对 100 个 数据 点 
的 光 变 曲线 建 模 ， 在 128GB 的 内 存 下 搭载 主 频 为 2.10GHz 的 Intel(R) Xeon(R) Gold 5218R 的 CPU， 
采用 该 模型 对 光 变 曲线 进行 建 模 , 需要 的 运行 时 间 约 为 Sms, 采用 PHOEBE 在 同样 的 运行 条 件 下 ， 
对 100 个 数据 点 的 拟 合 所 需要 的 运行 时 间 约 为 4s。 但 是 ， 如 果 采 用 MCMC+PHOEBE 对 一 条 光 变 
曲线 进行 参数 拟 合 ， 对 于 100 个 数据 点 ， 在 4 个 自由 参数 ，MCMC 的 walker 为 80，step 设 置 为 2000 
时 ， 则 需要 的 运行 时 间 约 10 hp 以 上 。 


4 总 结 和 展望 


本 文 基于 神经 网 络 构 建 了 半 接 双星 快速 测 光 解 轨 的 模型 ， 其 中 包含 主星 充满 和 次 星 充满 两 
测 光 解 轨 模型 。 模 型 根据 光 变 曲线 和 主星 温度 对 半 接 双星 的 光 变 曲线 进行 建 模 ， 从 而 获得 半 
接 双星 的 轨道 倾角 、 相 对 半径 、 质 量 比 和 温度 比 4 个 基本 参数 。 结 果 表 明 ， 针 对 模拟 数据 ， 主 星 
充满 测 光 解 轨 模型 对 轨道 倾角 、 相 对 半径 、 质 量 比 和 温度 比 的 残 差 的 标准 差分 别 为 0.238, 0.003, 
0.010 和 0.024， 拟 合 度 可 达 0.99; 次 星 充满 测 光 解 轨 模 型 对 轨道 倾角 、 相 对 半径 、 质 量 比 和 温度 
比 的 残 差 的 标准 差分 别 为 0.291，0.004，0.013 和 0.016， 并 且 拟 合 度 可 达 0.99。 模 型 对 Kepler 的 半 
接 双 星 解 轨 的 结果 表明 ， 对 于 对 称 的 半 接 双星 的 光 变 曲 线 ， 该 神经 网 络 构建 了 半 接 双星 快速 测 
光 解 轨 的 模型 对 参数 的 估计 较为 准确 ， 并 且 拟 合 度 可 达 0.9 以 上 。 在 相对 误差 为 0.01 时 ， 模 型 对 
于 轨道 倾角 接近 90”、 温 度 比 约 为 0.6、 光 变 振幅 为 1.84mag 的 半 接 双星 的 轨道 倾角 、 相 对 半径 、 


个 


= 


质量 比 和 温度 比 的 测量 误差 分 别 为 1.251，0.004，0.008 和 0.003。 另 外 ， 模 型 能 够 较为 准确 地 对 
脉动 食 双 星 的 光 变 曲线 进行 建 模 ， 这 对 后 续 分 析 恒 星 脉 动机 制 有 着重 要 的 意义 ， 其 次 ， 该 模型 
也 为 我 们 发 现 一 些 具 有 特殊 机 制 的 食 双 星 提供 了 一 个 有 效 的 工具 。 

采用 该 测 光 解 轨 模型 大 大 减少 了 光 变 曲线 建 模 需 要 耗费 的 时 间 ， 这 对 从 海量 光 变 曲线 中 获 
得 半 接 双星 的 轨道 参数 解 有 重要 的 意义 ， 这 些 参数 的 获得 为 研究 密 近 双星 系统 的 质量 和 角 动 量 
转移 、 损 失 和 重新 分 布 等 提供 了 有 利 的 数据 基础 。 并 且 ， 基 于 光 变 曲线 的 建 模 也 是 研究 食 双星 
中 的 恒星 脉动 机 制 的 基础 ， 通 过 光 变 曲线 建 模 消除 双星 掩 食 造 成 的 光度 变化 对 研究 恒星 的 脉动 
机 制 有 重要 的 意义 ， 这 为 我 们 了 解 恒星 内 部 提供 了 有 利 的 工具 。 同 时 ， 本 文 提 出 的 模型 针对 各 
个 望远镜 的 测 光 数据 均 可 进行 迁移 学 习 ， 例如， 如 果 针 对 TESS 的 数据 进行 分 析 ， 只 需 将 数据 和 
换 成 TESS 的 观测 波段 的 数据 集 进 行 训练 。 针 对 海量 的 光 变 曲线 ， 该 半 接 双星 的 光 变 曲线 建 模 工 
有 具 可 以 对 大 量 的 光 变 曲线 进行 快速 建 模 ， 这 对 开展 半 接 双星 的 性 质 统 计 分 析 以 及 快速 搜寻 半 接 
双星 中 具有 特殊 机 制 的 样本 有 重要 意义 。 
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Abstract: Semi-detached binaries are significant targets for the study of the formation and 
evolution of interacting binaries. Rapid modeling tool is highly required to derive the parameters 
with large amount of stars to be observed by many recent time-domain photometric surveys. In 
this work, based on a neural network, a light curve modeling of semi-detached binaries is 
proposed, which can derive orbital inclination(incl), relative radius(R/a), the mass ratio(q), and 
temperature ratio(T;/T;) fast via the observational light curve and known effective temperature of 
the primary star. The results of Kepler's light curve modeling show that the model can accurately 
fit the light curves of pulsating eclipsing binaries (the fitting degree can reach more than 0.9). For 
a target whose relative measurement error, orbital inclination, the amplitude of light curve, and 
temperature ratio are 0.01, ~90°, 1.84 mag, and 0.6, the measurement errors are 1.251, 0.004, 
0.008 and 0.003 for incl, R/a, q, and T;/T;, respectively. In addition, as an application, the 
proposed model in this work can be deployed on other photometric data by simply replacing the 
train data, which provides an effective tool to obtain a large number of parameters of 
semi-detached binaries and fast search for candidates of abnormal binaries. 
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